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Abstract 

We study the global dynamical behavior of spatially homogeneous solutions of the Einstein 
^ ' equations in Bianchi type I symmetry, where we use non-tilted elastic matter as an anisotropic 

CO I matter model that naturally generalizes perfect fluids. Based on our dynamical systems formu- 

. lation of the equations we are able to prove that (i) toward the future all solutions isotropize; 

00 ' (ii) toward the initial singularity all solutions display oscillatory behavior; solutions do not 

, converge to Kasner solutions but oscillate between different Kasner states. This behavior is 

• associated with energy condition violation as the singularity is approached. 

o 

'. 1 Introduction 



X 



Understanding the dynamics of cosmological models is one of the main goals of theoretical cos- 
mology. It is generally conceded that the analysis of generic cosmological models (which are 
space-times without symmetries) is exceedingly difficult; in fact, at present a mathematically rig- 
orous treatment seems to be out of reach. Nevertheless, heuristic and numerical studies that have 
been performed over the years have resulted in the formulation of a number of conjectures on the 
asymptotic dynamics of generic cosmologies, see |6j and references therein. In particular, it is 
conjectured that the generic singularity is intimately connected with (in fact, actually built on) 
the dynamics of spatially homogeneous cosmologies. 

Spatially homogeneous cosmological models have been analyzed intensively over the years, so that 
both the asymptotics toward the initial singularity and the asymptotics in the regime of infinite 
expansion (with infinitely diluted matter) are well understood; see [13] for a review. Most of 
the results concern solutions of the Einstein equations coupled to a perfect fluid which is usually 
assumed to obey a linear equation of state. However, it is unclear in general, how robust these 
results are under a change of the matter model. For example, it was shown in [5J that the structure 
of the initial singularity for coUisionless matter is considerably different from that of a perfect fluid 
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already for models of Bianchi type I. It is assumed that such different behavior stems from the 
anisotropic character of the stress-energy tensor |TT] . 

In this paper we investigate the global dynamics of spatially homogeneous solutions of Bianchi 
type I with anisotropic elastic matter. Elastic matter is described by the general relativistic 
theory of elasticity that was formulated by Carter/Quintana in [3] and further elaborated by 
Kijowski/Magh [7], Beig/Schmidt |2j and Karlovini/Samuelsson [8]; see also (TOlIIllIll]- For very 
recent work on the static Einstein-elastic matter equations see [1]. Our choice of matter model 
is motivated by the fact that elasticity theory offers a natural way of generalizing perfect fluids 
to a class of anisotropic phenomenological matter models without the need to resort to ad hoc 
assumptions on the expression of the anisotropic pressures. We consider elastic matter with a 
simple constitutive equation (Lagrangian) that leads to a stress-energy tensor of the form 

T^. = T^;;'^ + bU,,, (1) 

where T^""^ is the stress-energy tensor of a perfect fluid with linear equation of state and 6 is a 
constant which modulates the contribution of the anisotropic stress tensor n^,y. 

The main results of the paper are the following: Toward the future we observe isotropization of 
models. All solutions resemble infinitely diluted perfect fluid solutions in the asymptotic regime, 
which is in accord with physical intuition. Toward the initial singularity the behavior of Bianchi 
type I models with elastic matter is significantly different from the behavior of perfect fluid models. 
We prove that the behavior toward the singularity is oscillatory. In particular, there does not exist 
any solution that converges to a Kasner solution (vacuum solution). In the LRS case solutions 
oscillate between two different Kasner states (the Taub solution and the non-flat LRS solution); 
in the general case, the solution undergoes a (probably infinite) sequence of Kasner states (vac- 
uum states) as the singularity is approached. This behavior is intimately connected with energy 
condition violation. 

The paper is organized as follows. In Section[2]we briefly outline the derivation of the stress-energy 
tensor for elastic matter and introduce the class of diagonal Bianchi type I solutions of the Einstein 
equations which will be the subject of our analysis. In Section [3] we reformulate the Einstein 
equations with elastic matter for diagonal Bianchi type I solutions as a reduced dimensionless 
dynamical system on a compact state space. Section |4] contains the basic qualitative results on 
the global dynamics of solutions. In Section [5] we specialize to the LRS (locally rotationally 
symmetric) case; the reduced number of degrees of freedom permits a comprehensive analysis of 
the past attractor. In Section [6] we discuss another (non-generic) subclass of solutions ("anti-LRS 
solutions" ) whose behavior resembles that of LRS solutions. Finally, the most difhcult problem is 
addressed in Section [7l We present a detailed analysis of the past attractor of the full system and 
the associated past asymptotic behavior of generic solutions. 

2 Bianchi type I elastic spacetimes 

We begin this section with an introduction to the general relativistic theory of elasticity. However, 
since the only object of the theory used in this paper is the stress-energy tensor T^j, for elastic 
matter, we shall restrict ourselves to a discussion of the basic concepts needed for the derivation 
of T^jy. Comprehensive presentations of relativistic elasticity can be found in the references listed 
at the end of the paper. Most of the conventions we adopt, in particular those for the definition of 
the shear scalar and the elastic equation of state, are taken from [8]. (An option for the reader who 
is not interested in the derivation of T^j/, is to simply take as a given anisotropic stress-energy 
tensor and to proceed to the discussion of the Einstein equations in Bianchi type I, see ([8|) together 
with Uni) and dni).) 
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Let {M,g) denote the space-time, i.e., a four-dimensional manifold M with Lorentzian metric g 
that is time-orientable; local coordinates on M are given by x'', /j, = 0, . . . , 3. The material space 
(or body manifold) (iV, 7) is a three-dimensional Riemannian manifold; local coordinates are X^, 
A = 1,2, 3. Points in the material space identify particles of the material (in the continuum limit), 
where 7 measures the distance between the particles in the natural (unstrained) state of the matter. 
In the applications, 7 is usually chosen to be flat. The configuration function ip is defined to be a 
(smooth) map 

such that the kernel of the deformation gradient Tip : TM TN is generated by a (future-directed 
unit) timelike vector field u, i.e., kerT?/; = {u) or u^d^il)^ = 0. The vector field u is the matter 
four- velocity; by construction, ip~^{p) (i.e., the world-line of the particle p G N) is an integral 
curve. 

We define two metrics on the orthogonal complement (u)^ of u in TM (which coincides with TE if 
u is hypersurface orthogonal and thus generates a family of spacelike hypersurfaces E C M) . The 
Riemannian metric induced by g we denote by g: 

g — 9^1^ ~^ 'a^Ui/ . 

The pull-back of the material metric by the map -0, i.e., '0*(7), is called the relativistic strain 
tensor h: 

hfj^i, = df,il)^di,%l)^ JAB ; 

since hf^^u^ — 0, it is a metric in (u)^; since Cuh^y = 0, it is constant along the matter flow. The 
metric h^^ on (u)^ is Riemannian, hence /i^^, has three positive eigenvalues hi, h2, h^. 

The material is unstrained at x iff g^i,(x) = h^j^^(x). The scalar quantity 

is the particle density of the material. This interpretation is justified by virtue of the continuity 
equation 

(nw^) = . 

A specific choice of elastic material is made by postulating a constitutive equation, i.e., the func- 
tional dependence of the (rest frame) energy density p of the material on the configuration map, 
the deformation gradient and the space-time metric. An important class of materials is the one for 
which this functional dependence enters only through the principal invariants of the strain tensor. 
In this case we have 

p = p{qi,q2,qz), (2) 

where 

qi^tvh, q2 = tr[h^), 53 = tr (/i^) ; 

since = {qf ~ 3qiq2 + 2q3)/6, one of the invariants qi can be replaced by the particle density 
n. The materials described by ([2]) generalize the class of isotropic, homogeneous, hyperelastic 
materials from the classical theory of elasticity, see [9]. In many astrophysical applications (e.g., 
for the description of the solid crust of neutron stars), the effect of very large strains can be modeled 
by an elastic material in the quasi Hookean approximation [S]. This corresponds to a constitutive 
equation p that (i) depends on n; (ii) depends linearly on a quadratic invariant of the strain; (iii) 
has an absolute minimum at zero strain. Following |8] we choose the quadratic strain invariant to 
be the shear scalar, which is given by 

= 4 - - 24] , (3a) 
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or, in terms of the eigenvalues hi, /12, /13, 
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(3b) 



Evidently, is non-negative, and — (no shear) iff /i^^ cx g^^ (or equivalently, hi = h2 = h^). 
In this paper we shall consider a constitutive equation of the form 

p = /5(n)+/^(n)s^ (4) 

where p{n) is the unsheared energy density and jj.{n) the modulus of rigidity. The stress-energy 
tensor associated with these materials is obtained as the variation with respect to the space-time 
metric of the matter action Sm = — J \/W\ P- The result is given in [51 Sec. 6] and reads 



1 il 

where T^^, = P 9,j.u + 



- (tr(/i3) _ iiihf) g^^ + {trhfh^^ - [h^)^^ 



Here p is the isotropic (component of the) pressure, which is given by 



p = p{n) + v{n)s'^ 



where p = n — — — 
an \ n 



dp, 

n- /i 

an 



(5a) 

(5b) 

(6) 



The principal pressures Pi (which are the [non-zero] eigenvalues of T'^) are thus of the form pi = 
p + Spi] for an unstrained configuration, p^ = p, i = 1,2,3. For = (or = 0), the elastic 
material reduces to a perfect fluid with stress-energy tensor T^^ — pu^u^ +P9nu, energy density 
p = p and pressure p = p. 

It remains to specify the functions p and p, in the constitutive equation ([4]). We postulate a linear 
equation of state between the unsheared pressure p and the unsheared energy density p. 



p = ap 



(ae [-1,1]) 



and a linear equation of state between the modulus of rigidity p and the unsheared pressure p, 



p = bp {ab > 0) . 

By ^ this is equivalent to setting 



p^pon'^+'^, p^poahn°-+'^ (|a| < 1, a6 > 0) 

for some constant po > 0. Accordingly, 

p = pon^+^l + afos^) , p = ap. (7) 

Since for an unstrained material p — p and pi — p = p hold, i = 1,2,3, the bound |a| < 1 
ensures that the dominant energy condition < p is satisfied for an unstrained configuration. 
Furthermore, ab > guarantees that the energy density is positive for all values of the shear scalar 
and has a minimum at zero shear. When 6 = 0, the modulus of rigidity p vanishes and the 
elastic matter reduces to a perfect fluid with a linear equation of state p = ap; the condition |a| < 1 
ensures that the dominant energy condition \p\ < p is satisfied for this perfect fluid. When a = 
(so that p = 0), the choice of b is irrelevant, since ab = 0; this is clear because shear cannot occur 



4 



for dust. Henceforth, unless stated otherwise, by elastic matter we will always mean matter with 
constitutive equation (O, where a G [—1, 1] and ab > 0. 

Consider now a homogeneous space-time {M,g) of Bianchi type I, i.e., 

gfj^^dx'^dx" -dt^ + gij{t)dx^dx^ , (8) 

where gij{t), i,j = 1,2,3, is a family of Ricmannian metrics that is induced on the spatially ho- 
mogeneous hypersurfaces t — const. Let ^j;^{t,x^) : {M,g) {N,Sab) be a material configuration 
and df^Tp^ the corresponding deformation gradient. Compatibility with Bianchi type I symmetry 
forces the deformation gradient 9^-0"^ and thus the matter four- velocity to be functions of t only. 
We assume non-tilted matter: is orthogonal to the surfaces t = constant, i.e., = dt- This 
implies that = u^^dp^ip^ = dttp"^ and thus dtd^^^ — 0; hence dfj^ij)^ is constant with 9oV'^ = 0. 
For the strain tensor we find 

^00 = = 0, hi-j = SAsdiTp^djip^ = const ; 

since hj = g^^hjk, the matrix /i'^ is time-dependent as are its eigenvalues /ii, /i2, h^- It follows 
from ([SJ that 

Tm = P, Tfik = jk = 0, fij = (9) 
where is given in terms of hij via (|5b[) . 

The Einstein equations, in units c = 1 = SttG, decompose into the momentum constraint jk — 0, 
which is automatically satisfied by ([9]), the Hamiltonian constraint 



{trky - k',k^i - 2p = , (10a) 



and the evolution equations 



dtg.j = -2hj dtk\^ = (trfc)fc'^. - r^. -f -5'^ {T\ ~ p). (10b) 



Here, kij is the second fundamental form of the surfaces t — constant; Latin indexes are raised and 
lowered with gij . 

The Cauchy data associated with this initial value problem is given by (7y (0), fc*j(0); in addition 
we prescribe hij = hij{0) — const. Without loss of generality we can assume that gij{0) and fc*j(0) 
are diagonal (by choosing coordinates adapted to an orthogonal basis of eigenvectors of k'^j{0)). 
Furthermore we impose the condition that hij is diagonal; in particular, by rescaling the spatial 
coordinates, we can assume hij — Sij. This assumption is consistent with the evolution equations: 
Since the off-diagonal elements of the tensor T'^ form an homogeneous polynomial in h'^j = g'^^hjk, 
i j, it follows from the evolution equations (jlObp that {gij, k'^j, h^ j) remain diagonal for all times. 
Henceforth, we will restrict our attention to this special class of solutions of the equations pO|) . 
which are referred to as diagonal models. 

From /i* ■ = g^^hkj = diag((7^^, g^^, g'^'^) ~ diag(ft,i, /i2, ^3) we conclude that 
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1- 
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1- 

9 



- 6 



(11) 



cf. pb|) . which can be inserted into ([7]), i.e., 

P = Po(5"5'V')^''+'^/'(l + a^^s2), (|a|<l,a5>0), (12a) 
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to yield p as a function oi g^^, g^^, g^^. Moreover, from (jSbp we find 



T\^P+lf'(^^-^ + ^^^] ' (12b) 

1 . fg^^ 



g33 


5" 


5^' 




^ g22 




.9" 


g22 


^33 


g22 


%33 


^22 


g22 


^33 




g33 


^9" 


^33 



7^2 =P+ 0^(^^-32+33-722)' (12c) 
3 1 ^5'^ 



^3=P+«A(-22-33+ — -33)' (12d) 



where p — ap and /i = poa&(5^^3^^3'^^)^''^^-'^^ and are thus functions of g^^, g^^, g'^^. In the 
following we analyze the equations pO)) with anisotropic stress-energy tensor (fT2|) . In the diagonal 
case we consider, the unknowns are the six variables (g",/c*j) (no summation over i); the Cauchy 
data is {g" (0) , k\{0)) . 



3 Dynamical system formulation 

In order to formulate Einstein equations with elastic matter in Bianchi type I as a regular dynamical 
system we introduce alternative variables and matter quantities. Let 

i/ = -^, (13a) 

^^ = - 1 (no sum) ( El + E2 + S3 = ) . (13b) 

The Hubble scalar H carries dimension, while the shear variables Ei, E2, E3 are dimensionless. 
Evidently, the transformation between the variables /c^2: ^^3) and {H, Ei, E2, E3), where Si + 
E2 + E3 = 0, is one-to-one. 

In analogy to the new "momentum variables" we introduce new "configuration variables" . Let 
(ijk) be a cyclic permutation of (123). We define 

G = det5-i=giW', (14a) 

(eyfc = +1) \^ —- — — = 1 • (14b) 



5 ^ 1 - yi 1 - 2/2 1 - ys 

The variable G is dimensional, the variables yi, j/2, 1)3 are dimensionless; by construction, we 
have < j/i < 1 for all i. The transformation of variables (g^^, 5^^, g'^^) 1-^ {G,yi,y2,y3), where 
(yi, 2/2, 2/3) are subject to the constraint (yiy22/3)/[(l - 2/i)(l - y2)(l - Us)] = 1, is invertible, since 

{gn'-G^-y^ ie.,,^+l). (15) 

As a next step we normalize the matter quantities; we replace {p,T\,T\,T\) by (fi, tui, W2, ^3) 
which we define as 

T* p 
Wi = — - (no sum) , Q — 



P ' ' SH' 

It is customary to also introduce w by 



w ^ - y w, - ' — 



1 ^ 



3^ 3 p 
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Since p = ap for the elastic materials under consideration, see Section [2j we obtain 



Expressed in terms of the new variables the (dimensionless) shear scalar is given by 



1 
12 



3 

E 



1-% 
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(16) 



and the quantities Wi become 



Wi = a + 



ab \ yj 



Vi 



6 1 + a6 

where we have used the elastic equations of state of Section [5] 
Finally, we introduce a dimensionless time variable t defined through 



-1), 



(17) 



(18) 



and we henceforth adopt the convention that a prime denotes differentiation with respect to t. 
Written in the new dynamical variables the Einstein equations split into the dimensional equations 



H' 



-3H 



1 



G' 



-6G 



and a reduced set of dimensionless equations: 

'1 



-3n 



-(1 - - {w^ - w) 



y[ = -2yj(l - Vi) [Sj - Sfe] {eijk = +1) 
At the same time, the Hamiltonian constraint pOap reads 



(j = 1,2,3) 
(z = 1,2,3) 



1 



where S 



(19) 

(20a) 
(20b) 

(20c) 



We have thus arrived at the desired dynamical systems formulation of our problem: The dynamical 
system ([^0)1 describes the dynamics of Bianchi type I elastic spacetimes, where our choice of 
equations of state enters through the functions Wi(jji, 1/2,1/3), i = 1,2,3. Once the system ((20)) 
has been solved, the decoupled dimensional equations p9)l can be integrated and the standard 
variables, in particular the spatial metric, can be reconstructed. 

In addition to the dynamical system (j20p . the following auxiliary equation for f2 will prove to be 
useful: 



n' = n 



3(l-w)I]2-^u.feS, 



(21) 



In the remainder of this section we analyze in detail the state space X of the dynamical system ((20|l . 
The state space X is four-dimensional; it is defined as the Cartesian product of two two-dimensional 
sets, 

X ^Sxy, (22) 
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Si 




S2 2/2 
(a) The Kasner disc S (b) The space y 



Figure 1: The four-dimensional state space X is the Cartesian product of the Kasner disc S and the 
surface 3^. The space y is most conveniently represented as (the interior of) a hexagon; the center 
of the hexagon is the point (yi, 1/2, 2/3) = (1/2, 1/2, 1/2); for the edges the values of {yi, 1/2,2/3) are 
given in the figure, where the arrows denote the directions of increasing values (from to 1). 



where 

3 

r = {(Ei,S2,S3)| ^S»-0 A E2<i}^ (23a) 

i=l 

3 

:^' = {(2/i,2/2,2/3)|0<zj. <1V* A l[-y^ = l}. (23b) 

,=1 ^ 

The conditions on the variables are a direct consequence of the definitions (jl3bp and (jl4bp and the 
constraint (|20cp . 

The set S is the Kasner disc; it is usually depicted in a projection onto the plane with conormal 
(1, 1, 1), see Figured] The boundary of S is the Kasner circle KC = dS = {S^ = 1}. The Kasner 
circle contains six special points, which are referred to as LRS points: The three Taub points Ti, 
T2, T3 given by (S]i,E2,S3) — (2, —1,-1) and permutations, and the three non-flat LRS points 
Qi, Q2, Q3 given by (Si, E2, S3) = (—2, 1, 1) and permutations. The six sectors ofdS are denoted 
by permutations of the triple (123); by definition, < Ej < holds in sector (ijk). 

The set y is given as a two-dimensional surface in the interior of the unit cube [0, 1]^. Its boundary 
dy is the union of those six edges of the cube that do not contain the vertices (0, 0, 0) or (1, 1, 1). 
The projection of dy onto an affine plane with conormal (1, 1, 1) is a hexagon, the surface y itself 
its interior. The center of the hexagon represents the point (1/1,2/2, 2/3) = (1/2, 1/2, 1/2). In analogy 
to the Kasner circle, the six edges of the hexagon dy can be regarded as six sectors, where sector 
[ijk] is characterized by the inequality = yi < yj < yk = the six vertices of dy separate the 
sectors from each other: For % we have (y^, yj,yk) — (1, 0, 0), for Qi we have (y^, yj, y^) — (0, 1, 1); 
see Figure HJ 

A priori, by (|17[) . the quantities Wi are given as smooth functions of (yi,y2,y3) G y only. It is 
a crucial fact, however, that these functions admit a continuous extension to 3^, when we assume 
that ab ^ 0. In this case it is straightforward to prove that on sector [ijk] of dy, {wi,Wj,Wk) is 
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given by 

= a + 2ey7c(l - , = a - 2eyfc , Wk ^ a + 2eijkyj . (24) 

It follows that the r.h. side of the dynamical system ([20|) possesses an extension to the boundary of 
the state space, whereby we obtain a dynamical system on a compact state space X. In particular, 
the analysis of the flow on the boundary dX, which is 

dx = (ds xy)u(sx dy) , 

will turn out to be essential for an understanding of the global dynamics of the dynamical system. 

We conclude this section with some remarks on energy conditions. The dominant energy condition 
is expressed in the new matter variables as 

\w,\^l Vi = 1,2,3; (25) 

the weak energy condition reads 

-Is^Wi Vi = l,2,3; (26) 

the strong energy condition is satisfied if (I26p holds and a>— 1/3. It is a simple consequence 
of (jl7p that the dominant (and thus the weak) energy condition is satisfied for perfect fluids, i.e., 
for 6 = 0. However, for elastic matter, when \b\ > 0, the dominant energy condition is violated for 
(2/1,2/2,2/3) G y sufficiently close to the boundary dy. In fact, by ([24]) . on each sector [ijk] of dy 
there is at least one quantity Wi {i = 1,2,3) such that \wi\ > 1, and by continuity this inequality 
must hold in a neighborhood of that sector. (Note in this context that ([24ll is independent of the 
value oi b ^ 0.) On the other hand, provided that \a\ < 1, the dominant energy condition holds 
in some region of the interior of the state space, since Wi{l/2, 1/2, 1/2) = a for all i and thus by 
continuity {wil < 1 in a neighborhood of this point. 



4 Global dynamics 

The dynamical system (j20[) possesses one single equilibrium point in the state space X, which we 
call F. This fixed point is given by 

F : Si = E2 = S3 - , yi = y2 - 2/3 - ^ ; (27a) 

an alternative characterization is 

F : Si = S2 = S3 = , wi = W2 = = w = a , (27b) 

which is a direct consequence of (jl7p by taking into account the positivity of the variables yi. Since 
the principal pressures coincide, Pi = P2 = Ps = P, the fixed point F represents the flat isotropic 
FRW perfect fluid solution associated with the equation of state p = ap. 

We now consider the function 

Af = (1 - S^)-! (1 + ab s^) , (28) 

where is given by ((T6| and thus is a function of (2/1,2/2,2/3)- Recall that we suppose ab > 0. 
Accordingly, the function M is positive on the state space X] in fact, min;f M = 1 and the 
minimum M = 1 is attained at the fixed point F only. 
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A lengthy but straightforward computation, where we use the dynamical system ([20|) and the 
functions Wi{yi,y2,y^), see ([T7|. leads to 



on the subset = of the state space. 

Equation (|29ap suggests that a = 1 is a special case for our considerations, since then the r.h. side 
vanishes. We thus distinguish two cases: a < 1 and a—1. 

The case a < 1 

When we consider an elastic equation of state with a < 1 (and a6 > 0), Equation (|29p implies 
that M' < when 7^ and M"'\^2^q < except at the point F; therefore, M is a strictly 
monotonically decreasing function on X\F. The existence of a monotone function allows us to 
prove two central theorems. 

Theorem 1 (Future asymptotics). All orbits in the state space X converge to the fixed point F 
when T — !■ +00. 

Interpretation of the theorem. The fixed point F corresponds to a FRW perfect fluid solution 
associated with the equation of state p = ap; the theorem thus states that each Bianchi type I 
model with elastic matter obeying an equation of state with a € [—1,1) and ab > isotropizes 
toward the future and behaves like an (infinitely diluted) isotropic perfect fluid solution in the 
asymptotic regime. 

Proof. The function M is strictly monotonically decreasing along every orbit in the invariant set 
,Y\F. The monotonicity principle [U [13] implies that the w-limit of every orbit must be contained 
on the boundary, which is dX U F. Since M — +00 on dX, the boundary dX is excluded, which 
leaves the fixed point F as the only possible cj-limit. □ 

Theorem 2 (Past asymptotics). The a-limit of every orbit in X\F is a subset of the boundary 
dX of the state space. 

Proof. The monotonicity principle implies that the a-limit must be contained on dX U F. The 
point F is excluded, however, since M = minx Af = 1 at F. □ 

Theorem [2] describes the behavior of Bianchi type I models toward the initial singularity (which we 
choose to be t = 0). To see this we simply note that the inequalities —SH < H' < —3/2(1 + a)H 
follow from (fT9|) . so that p8|) can be integrated to yield a positive function ^(t) that satisfies t — > 
as T — > — 00. (In this context it is necessary to assume a > — 1; the case a = — I requires a different 
argument that involves more detailed information on the a-limits of orbits in X, see Section [T]) 

Since the boundary dX contains the Kasner circle(s), see Section [71 the theorem suggests that 
the Kasner solutions will play an essential role in the past asymptotic dynamics of Bianchi type 
I solutions with clastic matter. However, in Section [7| we will prove that there does not exist any 
solution that converges to a Kasner solution as t ^ 0; instead we observe oscillatory behavior 
toward the singularity. In order to derive these results concerning the past asymptotic dynamics 
of solutions we must analyze the flow of the dynamical system on the boundary dX , which will be 
done in Section [7] 



M' = -3(1 - a)j:^M . 



(29a) 



The computation of higher derivatives reveals that 




(29b) 



k 
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The case a = 1 

When a — 1, the function M is constant along the orbits of the dynamical system, i.e., for every 
K 9 m > 1, the hypersurface 

Mm = {(Si, S2, S3, yi, 2/2, yg) e A- 1 M = m > 1} 

is an invariant subset in A". In other words, M acts as a conserved "energy"; solutions with the 
"energy" M — m are contained on When to = 1, we have Mi — {F}. 

It is not difficult to show from (fT6|) that each surface = const is a (topological) sphere in the 
unit cube (0, 1)^, which is centered at (2/1,2/2, 2/3) = (1/2, 1/2, 1/2). Accordingly, — const defines 
a closed curve in y. Consequently, when we rewrite (|28p in the form 

MY? +ahs^ = M ~l (M = to > 1) , (30) 

we conclude that Mm represents a topological hypersphere in X whose center is the fixed point 
F. In particular. Mm does not intersect dX. We have proved the following result. 

Theorem 3. For the uj-limit set a;(7) of an orbit 7 in X\F we have: F ^ w(7) and u}{'^)C]dX = 0. 
An identical statement holds for the a-limit set ofj. 

Interpretation of the theorem. If the equation of state of the elastic matter is such that a = 1, 
then the associated Bianchi type 1 solutions of the Einstein equations do not isotropize toward the 
future (and neither toward the singularity). Furthermore, the solution cannot be approximated by 
Kasner solutions at any time and neither asymptotically (since Kasner solutions are represented 
by points on dX, see Section [T]). 

From (|30p it follows that the maximum shear of a solution with "energy" Af = to is given by 
— (to — l)/{ab), while remains bounded by < 1 — 1/to, so that f7 > 1/to. 

5 LRS solutions 

In this section we consider a special case of Bianchi type I models: Locally rotationally symmetric 
(LRS) models. The additional symmetry that is imposed reduces the number of degrees of freedom, 
which facilitates the analysis of the (past asymptotic) dynamics of solutions in both cases a < 1 
and a = 1. 

A Bianchi type I solution of the Einstein equations is locally rotationally symmetric (LRS) if 
Ej = Efc and Wj = Wk for some pair j 7^ k. Let (ijk) be the completion of the pair (j, k) to 
a permutation of the triple (123). The first condition, i.e., Ej = E^, implies that 2/i = const 
via (|20bp . and thus g-'-' /g^^ = const according to ([l4b|; by a possible rescaling of the spatial 
coordinates we obtain LRS geometry, i.e., g^^ = g'^'^ . The second condition, i.e., Wj = Wk, states 
that the matter content is compatible with LRS symmetry; it guarantees that j/^ = 1/2 (so that 
gjj = gkk automatically). To see this we use (fTT]) and observe that Wj — Wk iff 

1 - 2/j _ Vj ^ ^-yk _ 2/fc ^ 2 1 - y» _ Vi \ 
Vj 1 - Vj Vk 1 - 2/fe \ Vi 1 - Vi ) 

multiplication with (2/j2//c)/[(l ~ 2/i)(l ~ Vk)], where we use the constraint (jl4bp on the variables 

{yi,yj,yk), yields 

( + i/fc \ A _ Vj Vk \ = _2 fi - y'' \ (i + _Mj__y^\ 
Vi-% i-yfe/V 1-yji-y/c/ \ i-yji-2/fc/\ i-yji-yfe/ 
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Figure 2: In the state space X = E x y there exist three invariant subsets that can be identified 
as LRS subsets. We depict the intersection of LRSi, LRS2, LRS3 with the factors S and y of the 
state space. 



and hence 



Wj — Wk 1 — 



Vk 



Q ^ Vj ^ I - Vk ^ Vi 



In the state space X, the conditions S^- = and Hi = \ wj = Wk) define three invariant subsets 
which we denote by LRSi, « = 1,2,3, see Figure [5J In the dynamical systems representation of 
Bianchi type I elastic cosmologies, LRS configurations are given by orbits on one of these invariant 
subsets. In the following we thus analyze the dynamics of LRS solutions by studying the flow of 
the dynamical system on the LRS subsets LRS;. 

Let again (ijk) be a cyclic permutation of (123). Consider the subset LRSi C X, which is given 
by the conditions 

yi = ^ ( <^ Wj = Wfe % = 1 - Vk) and = Sfe ( 4^ -2Sj = -2Sfe) . 

On LRSi, the shear scalar reads 

2 _ 1 (1 - 2y,)2 _ 1 (1 - 22/fe)2 



6 - Vj) 6 yk{l - Vk) 



and the matter quantities become 

ab 



1 - 2y, 



6 y,{l-y,) + f{l-2y,y 
The dynamical system ([20)1 reduces to 



Wk 



Wi = 3a — 2wj 



V'j = J'jO^i 



-3(1 



-(l-a)E, 



[wj — a) 



(31) 

(32a) 
(32b) 



where we have used that Q, = 1 — Y?, under the present assumptions. 
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The state space LRSi can be represented as (the interior of) the rectangle (—1, 1) x (0, 1) 9 (Ej, y^). 
Since wj extends smoothly to yj = and yj = 1, the dynamical system (|32p extends smoothly to 
the compact space [—1,1] x [0, 1]. The four sides of the rectangle are invariant subspaces; when we 
exclude the vertices from our considerations, we find: 

. S,=±l y;. = T6y,(l-t/,) ^0 

m yj = => E^- = -|r2[(l-a)Sj+2] < 
• yj ^ 1 => E^- = -|r2[(l-a)Ej -2] > 

The four vertices of the rectangle are fixed points. (The notation is chosen to agree with the 
conventions of Section [T] ) 



Fixed point 


(Ei, Ej, Efc) 


{yt,yj,yk) 


Fixed point represents 


Q(jifc) 
Q(kij) 


Q, = (-2,+l,+l) 
Q, = (-2,+l,+l) 
T, = (+2,-1,-1) 
T, = (+2,-1,-1) 


(1/2,0,1) 
(1/2,1,0) 
(1/2,0,1) 
(1/2,1,0) 


non-flat LRS Kasner solution 
non-flat LRS Kasner solution 
Taub solution (flat LRS Kasner) 
Taub solution (flat LRS Kasner) 



From this analysis it follows that the boundary of the state space LRSi forms a heteroclinic cycle: 



^ (kij) 



Q{kij) 



(33) 



In the interior of the state space LRSi there exists one single fixed point: The FRW perfect fluid 
fixed point F; recall that E^ = Ej = E^ = and yi — y.j = = 1/2 at F. 

To analyze the global dynamics on LRSi we distinguish the cases a = 1 and a 7^ 1 as in Section [H 
First, let a = 1. The considerations of Section 3] imply that there exists a family of invariant 
subsets {Mm \ m > 1}. For each m, the surface M-m is a three-dimensional hypersphere and 
its intersection with the two-dimensional surface LRSi yields a closed curve with center F. By 
construction, this closed curve is an orbit of the dynamical system on LRSi. More explicitly, we 
see that 

Af = (1 - + abs^) = (1 - Y^r' f 1 + ^ ^^77^) = - > 1 (34) 

defines a family of periodic orbits in LRSi which are centered at the fixed point F. The phase 
portrait of the dynamical system in the case a = 1 is represented in Figure [3l (Note that in the 
limiting fluid case, i.e., a = 1 with 6 = 0, the orbits are no longer periodic but straight lines 
Ej ~ const.) 

Second we consider elastic matter with the property a 7^ 1. The global dynamics on LRSi follows 
from Theorem [1] and Theorem [2] in conjunction with the analysis of the boundary of LRSi: 

Corollary 4. The uj-limit of every orbit on LRSi is the fixed point F, the a-limit is the heteroclinic 
cycle ([33|. 

Interpretation of the corollary. Each LRS Bianchi type I model with elastic matter obeying 
an equation of state with a G [—1,1) and a6 > isotropizes toward the future and behaves like an 
(infinitely diluted) isotropic perfect fluid solution in the asymptotic regime. Toward the singularity 
we observe oscillatory behavior between the non-flat LRS Kasner solution and the Taub solution. 
This is in stark contrast to the behavior of perfect fluid solutions, which converge to either the 
non-flat LRS Kasner solution or to the Taub solution as t ^ 0. 
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\jik) Q(jife) 



Figure 3: Phase portraits of LRS solutions for a ~ 1 (left) and a G [—1, 1) (right). The dominant 
energy condition is violated in the shadowed region (for a = — 1, this is everywhere except at 
yj = 1/2). It is assumed that ab < 3 and, for a E [—1,1), that ab > -^{1 — a)^ (oscillatory 
approach toward F). 



The past asymptotic dynamics of solutions as described by the corollary is intimately connected 
with the violation of energy conditions. As a matter of course, the general statement of Section [3] 
also applies in the LRS case: While the dominant energy condition (and thus the weak energy 
condition) is satisfied in a neighborhood of the fixed point F (where it is assumed that \a\ < 1), 
we observe energy condition violation in a neighborhood of the boundaries yj = and yj = 1. 

The discussion of the issue of energy condition violation is facilitated when we assume the upper 
bound ab <3. Under this condition it is straightforward to show that {wj — a) is a monotonically 
increasing function of yj £ [0, 1] with range [—1, 1]; accordingly, {wi — a) is decreasing with range 
[—2,2]. A particularly simple case, which displays all the relevant features of the general case 
ab < 3, is ab = 3/2, since (wi, Wj, Wk) become linear in yj, i.e., 

Wj ^ a + 2yj — 1 = Wk , Wi ^ a + 2 — Ayj . (35) 

A straightforward calculation yields the following results for this special case: 



The weak energy condition is satisfied iff 



a a + 3 

'r 4 



This interval collapses to the point yj — 1/2 when a = ~l. 
The dominant energy condition is satisfied iff 

a + 1 a + 3 



a a 

r'^2 



n 



4 ' 4 

see Figure |3l This domain reduces to the point yj ~ 1/2 when a ~ ±1. 



A general property that is worth observing is the presence of regions where the dominant energy 
condition is violated, while the weak condition is satisfied. 

For a given a £ (—1,1), the regions of energy condition violation become smaller when we let 
|6| — > 0. However, the general statement of Section [3] applies for all b, no matter how small: There 
exists a neighborhood of the boundaries yj = and yj = 1, where the energy conditions are 
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violated. To first order in 6, the region where the energy conditions hold is characterized by 

{l + a)ab a^b 1 f (l - a)ab a^b 
< 2/i < nun < 1 — — , 1 



6(1 -a)' 6(l + a)j-"^- \ 6(1 + a)' 6(1 - a) 
only when 6 = 0, i.e., in the fluid case, we obtain < j/j < 1. 

To conclude this section we study in detail the future asymptotics of solutions, which in the state 
space description of the dynamics corresponds to investigating the flow in the neighborhood of the 
fixed point F. 

Since the r.h. side of the dynamical system is smooth in a neighborhood of F we can perform a 
local dynamical systems analysis. The linearization of the dynamical system at F possesses the 
eigenvalues 



\i^\[-{l-a)-sj{l-aY~^abj , = ^-(1 - a) + ^(1 - a)^ - | 

The eigenvectors associated with Ai^2 are 



Via = i ab 
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(^l_a)T\lil-ar~fab 



(36) 



It is immediate that 



• the eigenvalue(s) are real (and negative), if ab ^ -^{1 — a) ; in this case F is a stable node; 

• the eigenvalues are complex (with negative real part), if ab > ]^(1 — a)^; in this case the 
fixed point F is a stable focus and the solutions' approach to F as r — > oo is oscillatory; see 
Figure [31 

The late time behavior of Bianchi type I models with clastic matter is thus characterized by 

• monotonic isotropization if ab ^ (3/32)(l — a)^; 

• oscillatory isotropization if ab > (3/32)(l — a)'^; both the amplitude of the oscillations and 
the frequency are decreasing as t oo. 



6 Anti-LRS solutions 

The LRS subsets of the Kasner disc S are defined by requiring that = Efe for some pair {j, k). 
Analogously, we define the three anti-LRS subsets by setting Y.j = — Sfc for some pair {j,k); 
consequently, = 0, where (ijk) denotes the completion of the pair (j, k) to a permutation of 
the triple (123). Let eijk = 1; the six anti-LRS points on the Kasner circle dS are given by the 
three points (I]i,Sj,Efe) = (— V3, 0, VS) G i^jk), which we denote by Pj, j — 1,2,3, and by the 
three points (E^, Sj, E^) — (\/3, 0, — \/3) G {kji), which we denote by Sj, j = 1, 2, 3. (The notation 
is chosen to complement the standard notation Qj, Tj for the LRS points.) Since the angular 
distance on the Kasner circle is a nai've measure for the difference between Kasner states, anti-LRS 
states are those which are maximally different from LRS states; hence the name. 

Let again {ijk) be a cyclic permutation of (123). We define the anti-LRS subset ^LRS^ in the state 
space X by 

"LRS, : = -Sfe A = yk ■ (37a) 
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A equivalent definition is 

'^LRS, : = A w^=w = a; 

1,2,3, so that there exist three anti-LRS subsets: ^LRS;, i 



(37b) 



here, i = 1,2,3, so that there exist three anti-LRS subsets: ^LRS;, i = 1,2,3. Since = 
and {uj — Uk)' = on '^LRSi, cf. (|20|) . these subsets are invariant subsets in X. Orbits on '^LRSi 
generate anti-LRS solutions of Bianchi type I; the metric of these models is characterized by the 
condition that gu be the geometric mean of gjj and gkk, i-e.. 



9ii - \J 9jj 9kk ■ 

Evidently, a solution is both LRS and anti-LRS iff it is isotropic. 
On the anti-LRS subset ''LRS^ the dynamical system (piU) reduces to 

y'j = 22/j(l - yj)T.j , 

^ = -3(1 - S2/3) i(l - a)E, - (w, - a) 



where we have used that i7 = 1 — S|/3. It is not difficult to show that Wj reads 

2a6(l-2y,)[l-y,(l-y,)] 



a + 



12yKl - y^)2 + abil - 2y,)2[l - 2y,(l - y,)] 



(38) 

(39a) 
(39b) 

(40) 



The (closure of the) state space '^LRSi can be represented as the rectangle [— v^, -v/S] x [0, 1] 3 
(T.j,yj). The four vertices are fixed points. 



Fixed point 


(Sj, Sj, Efe) 


{yi,y],yk) 


Interpretation 


Pr. = {Pi} X {TJ 
Pq. = {Pi} X {Q.} 
Sr = {SJ X {TJ 
Sq, ={S^} X {Q,} 


Pi = (0,W3,-V3) 
P, = (0,+V3,-\/3) 
S, = (0,-V3,+V3) 
S, = (0,-V3,+V3) 


% = {1,0,0} 

= (0,1,1) 

% = (1,0,0) 

= (0,1,1) 


anti-LRS Kasner solution 
anti-LRS Kasner solution 
anti-LRS Kasner solution 
anti-LRS Kasner solution 



The boundary of ^LRS^ forms a heteroclinic cycle, 

Sq, < Pq. 



Sr 



(41) 



and in the interior of the space there is the fixed point F, cf. ([55]) . 

The global dynamics of anti-LRS solutions is therefore reminiscent of the dynamics of LRS so- 
lutions. We merely state the results; the proofs are analogous to the proofs of Section [5] Each 
anti-LRS solution isotropizes toward the future. The eigenvalues of the linearization of the dy- 
namical system at the point F are the same as in the LRS case; therefore we distinguish two kinds 
of isotropization: Monotonic isotropization and oscillatory isotropization. Toward the singularity 
we observe oscillatory behavior between the two anti-LRS Kasner states; this regime is connected 
with energy condition violation. 



7 Past asymptotic states 

This section is devoted to investigating the dynamics of general diagonal Bianchi type I solutions 
toward the initial singularity. Theorem [2] states that (i) if the elastic matter satisfies a e [—1,1) 
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(which we will assume from now on) and (ii) if we exclude the isotropic FRW solution represented 
by F from our considerations, then the a-limit set of every orbit in the state space X is located on 
dX. Hence, in order to understand the structure of the a-limit set, it is necessary to study in detail 
the flow on the boundary. As a preparatory step, we discuss the network of fixed points (which 
includes the Kasner circles) that is present on dX . In Subsection 17. 11 we analyze in a step- by-step 
manner the flow on the invariant subsets of dX; in Subsection 17.21 these results are combined to 
identify certain special structures on dX , namely heteroclinic cycles and heteroclinic sequences. 
Finally, in Subsection 17.31 we condense the collected results into statements (such as Theorem [S]) 
and conjectures on the possible a-limit sets on dX. 

Since X — S x y, the boundary dX consists of two components, 

dx = {as xy)u{sx dy) , (42) 

the intersection of which is the set dS x dy. 

By construction, the boundary sets dS x y and S x dy are invariant under the flow of the (induced) 
dynamical system (j20p . There exists a number of equilibrium points of the dynamical system; these 
fixed points are located not in the interior of the boundary components (i.e., neither on dS x y nor 
on 17 X dy), but on the shared boundary dU x dy. In fact, on dE x dy there exists a connected 
network of one-parameter families of equilibrium points; note in particular that there do not exist 
isolated fixed points. 

• Kasner circles: There exist six families of fixed points that can be interpreted as Kasner 
circles (and are thus associated with the Kasner solutions). They arise at the six vertices 
{Ti,r2,T3, Qi, Q2, Q3} of dy, i.e., we have 

KCt, ■.= dS X {%} (fc = l,2,3), KCq^ :=ar X {QJ (i = 1,2,3). (43) 

• Taub lines and non-flat LRS lines: There exist two lines of fixed points associated with each 
edge (= sector) of dy. Consider the sector [ijk] = {(2/1,2/2, Vs) \ Vi ^ 0, < Vj < 1, Vk = 1}; 
then the Taub line and the non-flat LRS line associated with [ijk] are given by 

TL[,,.fc] {T,} X [ijk] , QL[,^.,] {Q,} x [ijk] , (44) 

respectively. (For later purposes, we also introduce the associated open sets TL(ijfc) and 
QL(ijfc)- While TLfijfe] = {T^} x [ijk], where [ijk] {0 = < yj < yu = 1}, T^l^(ijk) is 
given by TLjij^) = {Tj} x {ijk), where {ijk) = {0 = t/i < yj < yk — 1}, and analogously for 

• Since [ijk] connects Qi with 7^, the intersection of these families of fixed points (I43p and 
consists of 24 special points, which are {T^} x {7^}, {Qj} x {Tk}, {Tj} x {Qi}, {Qj} x {Qi}, 
where i ^ j ^ k. 

This fixed point structure on the boundary of the state space dX is depicted in Figure ID 

For the subsequent analysis of the flow on the components of dX we make the assumption a > — 1 . 
The case a = — 1 is degenerate in the sense that the linearized dynamical system at certain fixed 
points vanishes (i.e., all eigenvalues are zero). Since this leads to some technical difficulties that 
we do not want to discuss here, we henceforth assume a G (—1, 1) (and ab > 0). However, we note 
without proof that the key statements we derive in the following apply to the case a = — 1 as well. 
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Figure 4: A schematic depiction of the fixed points on the boundary dX. 



7.1 Invariant subsets 

The boundary component dS X y 

The set dS x 3^ is the subset of X characterized hy ft — (or, equivalently, = 1). The induced 
dynamical system on dU x y reads 

yl ^ -2eijkyi{l - yt)[T,j -Ek] (no summation) , (i = 1,2,3); (45) 

in particular, Sj = const for all i. Since = and — const, orbits on dS x y represent vacuum 
solutions (Kasner solutions); accordingly, the set dS x y can be called the vacuum subset. 

The space dS x y can be depicted as (set of points contained within or lying on) a torus whose 
cross section is the hexagon 3^. Since = const for all i, each cross section {(Si, S2, S3)} x y 
is an invariant subspace. Let (Ei, S3) be a element of sector (ijk) of dS, i.e., < < Efc. 
Then y[ cx Eijk, y'j oc {—eijk), y'k tijk, hence the variables yi are increasing or decreasing for all 
I. If tijk = 1, the a-limit of each orbit (in the interior of the space, i.e., in {(Si,E2,E3)} x 3^) is 
the point {(£1,^2,23)} x {7)} and the w-limit is the point {(Si, E2, S3)} x {Qj}; if Sijk = — 1, 
the roles of the points are interchanged. 

When (Si, S2, S3) is one the Taub points Tj, then orbits in the interior space {Tj} x 3^ emanate 
from a fixed point on the Taub line TL(j.ji) — {Tj} x (kji), where etjk = 1, and end at a fixed 
point on the Taub line TL(ijfe) = {T^} x (ijk). The result is converse, when (Si,S2,S3) is one 
of the non-flat LRS points Qj. In that case, orbits in the interior space {Qj} x y connect a fixed 
point on the non-fiat LRS line QLj-j^j,-) = {Qj} x (ijk), where e^fe = 1, with a fixed point on the 
Taub line Q'L^kji) = {Qj} (^i*)- 

The flow on the set dS x 3^ is depicted in Figure El 
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2/2 




S3 






tJ \ 


/ \Qi 


1 \ / \ / 



Figure 5: Flow on the boundary component dS x y. The cross sections of the "torus" dS x y 
are invariant subspaces, since (I]i,E2,S3) = const. The flow on the section {(Si, £2,23)} x y 
depends on the position (sector) of (Ei, £2, S3) on dS. Note that all fixed points on dS x y are 
hyperbolic or transversally hyperbolic. 



The boundary component S X dy 

Since dy consists of the six sectors [ijk], see Figure [1] the set S x dy can be viewed as the union 

Sxdy = \J(Sx [zjk]) = y Cyl[,^fcj . (46) 

ij k ij k 

When written out explicitly, we see that the set 

Cyl[,^fcj =r X {(yi, 2/2,^3) I = 0, < < 1, z/fe = 1} =^ X [ijk] (47) 

represents a cylinder, see Figure [6l The six cylinders Cyl^^j^ are aligned along the hexagon dy, 
where each vertex corresponds to the top/base of a cylinder: E x Qi {i = 1,2,3) and E x Tk 
{k = 1, 2, 3), respectively. 

The flow on the space Cyljjj;.] is given by the induced dynamical system 

= (48a) 
2eyfc y'j = -2eijkVj{l - yj){^k ~ ^i) (48b) 

yfc = 1 . (48c) 

Since (1 — a)|S]j| < 4 in the interior of Cyl^^j^ and in the interior of the top and the base of the 
cylinder (i.e., on i7 x [ijk]), the derivative of has a sign, so that Ej is a monotone function. 
The monotonicity principle thus implies that the a- and the w-limits of orbits must be located on 
the lateral surface of Cyljjjj,] (which includes the Kasner circles at the base/top). 



E', = -m 




a)Ei 


E^- = -sn 




a)Ej 




\>- 


a)Efe 
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KCr, 

Figure 6: The boundary component S x dy consists of six cylinders Cyljij^,] . The top of Cylyij^ = 
S X [ijk] is the Kasner disc S x Qi, the base of Cyl^.^j^ is the Kasner disc S x Tk. In this figure, 
Cyl[jjj.] (with djk = 1) is depicted together with the flow of the dynamical system on the lateral 
boundary. 



The flow on the lateral surface of Cyljjjj,] is simple. Since (Si, , Efc) = const, the equation 
y'j = —2eijkyj{l ~ y.j){Y,k ~ S^) contains the entire dynamical information. Suppose that e^fe = 1. 
Then, for Cyl^^jk] (and identically for Gy\\^kji]) ^nd y'^ ^ when ^ E^. Furthermore, y'^ — 
when Ei = E^, which is the case at the lines of fixed points TLjij^] and QLjjj-^.]. Therefore, the 
lateral surface consists of two domains where yj is increasing/decreasing and which are separated 
by the Taub line and the non-flat LRS line; see Figure [6l 

This analysis leaves the equilibrium points, i.e., the sets TL[,yfc], QlJ^^j^, KCg. , KCt^ as the only 
possible a- and w-limit sets of orbits in Cylj^^-j;.] . It thus merely remains to identify those equilibrium 
points that act as sources/sinks for interior orbits. 



(49) 



To this end we make use of the auxiliary equations l|2ip ; we find that 

n-^n' - 3(1 -a) + 2e,jk [(1 - - E,) + y,(Ej - E^)] 



when evaluated at a fixed point on the lateral surface. Consequently, for the equilibrium points 
we obtain 

(50a) 
(50b) 
(50c) 
(50d) 

L ^ J 

Suppose that = 1 [tijk = — 1]. Then each fixed point on TL^^jj.) acts as a source [sink] for 
one interior orbit. (Since the points on TL^^jj.) and QLj-j^-j,-) are not transversally hyperbolic, this 



TL[ij-fc] : 


n- 




3 

^ 2 




- eyfc)(3 - a) - 


(1-^ 


ijk){l + a) 


QL[ijfe] : 


n- 


-If!' 


3 

^ 2 


"(1- 


- eyfc)(3 - a) - 


(1 + ^ 


ijk){l + a) 


KCs, : 


n- 


-If!' 


= 2 


^(1 

2^ 


- a) + eijfe(Ej 


-Sfc) 


1 


KCr, : 


n- 


-If!' 


= 2 


^(1 


- a) + eijfe(Ej 







20 



KG 



Qi 




QL 



[ijk] 



KCr. 



3 



Figure 7: A schematic depiction of the flow on the cyhnder Cyl^^j^ = S x [ijk] {etjk = 1) 



terminology is to be understood in a broad sense. The precise statement reads: Each fixed point 
on TL(ijfc) acts as the a-hmit [w-hmit] for one interior orbit, while there do not exist any other 
orbits that converge to this point as t — > ±oo.) Analogously, the points on QL^j^-^-) act as sinks 
[sources] . 

The relations (|50cp and (jSOdp . in conjunction with the properties of the flow on the lateral boundary, 
lead to the following classiflcation of the flxed points on the Kasner circles: A flxed point on KCg^ 
is a (transversally hyperbolic) 

3 

source/sink ^ e„-fe(Ej - Sfe) ^ --(1 - a) and €ijk{^i - ^k) ^ . (51a) 

A fixed point on KCt^ is a (transversally hyperbolic) 

3 

source/sink ^ eyfe(Ei - S^) -(1 - a) and eyi;(Si - Efc) ^ . (51b) 

The remaining equilibrium points on the Kasner circles do not attract interior orbits as r — > ±oo. 
(Except for four special points, these fixed points are transversally hyperbolic saddles.) 

The results of our analysis of the flow on the set Cyl^^^j^ are summarized in Figure[71 Each interior 
orbit (i.e., each orbit in the interior of Cylj^j^.]) is a heteroclinic orbit; it has a source as its a-limit 
and a sink as its w-limit. 

We conclude our study of the flow on the set Cyl[jjfc] by noting that certain orbits can be given 
explicitly. First, there is the solution given by = Efc and yj = 1/2, which is LRS^-; this solution 
is a straight line that connects the fixed point T(ijk) — Tj x {(0, 1/2, 1)} on TL^^jk] with the point 
Q(ijfc) = Qi X {(0, 1/2, 1)} on QL[jj-j;,], see also ([55]) . Second, the orbits on the base/top of Cyl[jjj,], 
i.e., S X {Qi} and S x {Tk}, are of a simple geometric form. 

To see this, consider the top of Cyl^^j^, i.e., S x {Qi}- From the equations 



i(l - a)'Ej/k ± 2e^Jk 



(52) 



it is immediate that = is a solution of the system; clearly, this orbit is the intersection of 
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Figure 8: Flow on the top (left) and base (right) of Cyljjj^.] {eijk = 1). 



Cyl[jjj,] with the anti-LRS set '^LRSi. Furthermore, we find that 



ijk 



— constant (53) 



under the flow of the system, whereby we obtain an explicit representation of all orbits: Each orbit 
is a straight line. (Setting const = — 1 reproduces the anti-LRS orbit = 0.) The one-parameter 
family of straight lines described by Equation (|55|) possesses a common focal point, i.e., all lines 
intersect in the point (S^, Ej, E^) = j^eijk{0, —1, 1). Two members of this family of straight lines 
are tangential to the Kasner circle KCg. . The two associated points of contact lie on the straight 
line eijkCS^j — E^) = — 1(1 — a); these fixed points are not transversally hyperbolic; see (|51ap . For 
a depiction of the flow on S x {Qi} and the analogous flow on S x {7^}, see Figure[8l 

The three special orbits that exist on Cylj^j^.] are particularly relevant for further purposes. Let us 
thus recapitulate: There exists an anti-LRS orbit '^LRSi given by E^ = on the top of the cylinder, 
an LRSj orbit E^ — Y.^ in the middle of Cyljjj^,] , and an anti-LRS orbit '^LRSfc given by E^ = 
on the bottom. 



7.2 Structures on the boundary 

A crucial property of the equilibrium points on dX is that all these points are saddles. This is 
straightforward to prove: Consider first a fixed point on TL[ijfe] = {Tj}x[ijk]. As shown previously, 
such a point acts as a source [sink] within the subset Cyljjj^.] — S x [ijk] when eijk = 1 [e^fe = 
— 1]. However, within the (linearly independent) invariant subset {Tj} x 3^ it is a (transversally 
hyperbolic) sink [source]. Consequently, the point acts as a saddle. In particular, orbits in X 
cannot converge to TLj^jfej as r — > ±oo. The case QLj^j-^.] is analogous. 

Consider now the Kasner circle KCq. . This circle is part of the boundary of two independent 
invariant subsets: Cyljjjj,] — S x [ijk] and Cyljj^j] = S x [ikj]; without loss of generality we 
assume e^fe = 1. From (|5T|) we conclude that a point on 

3 

KCq, C Cyl[jjfe] is a source/sink <J4> (E^ - Ej.) ^ --(1 - a) , (E^ - E^,) ^ , (54a) 

3 

KCq, C Cyl[,fcj] is a source/sink ^ (E^ - Efe) ^ --(1 - a) , (E^ - E,) ^ . (54b) 

Since these conditions are mutually exclusive, every equilibrium point on KCg^ acts as a saddle. 
The analogous considerations apply to the case KCr^ ■ 
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Heteroclinic cycles 



From the analysis of the flow on the invariant subsets of dX we infer the existence of a multitude 
of heteroclinic cycles. On the one hand there are the three LRS cycles (|33|) . which arise as the 
intersection of dX with LRSi, i — 1,2,3; see Section [5l On the other hand there exist the three 
anti-LRS cycles (|4T|) . which are the intersection of dX with '^LRSi, i — 1, 2, 3, so that = 0; see 
Section [6l There are, however, infinitely many heteroclinic cycles associated with the condition 
Si = 0; this will be discussed next. 

On the boundary component dS x y consider, without loss of generality, the subset Si = 0. It 
consists of the two disconnected sets {Pi} x y and {Si} x 3^, where the point Pi = (0, \/3, — \/3) is 
the anti-LRS point of sector (312) and Si = (0, — -s/S, ^/3) the anti-LRS point of sector (213) of dS. 
The a-limit of each orbit in {Pi} x 3^ is the fixed point {Pi} x {Ti}, the w-hmit is {Pi} x {Qi}, 
see Figure [S] Note that these fixed points are not only connected with each other through orbits 
in {Pi} X 3^, but also via the sequences of boundary orbits 

{Pi} X {Ti} {Pi} X {Q,} ^ {Pi} X {%} ^ {Pi} X {Qi} 
and {Pi}x{ri}^{Pi}x{Q3}--{Pi}x{r2}^{Pi}x{Qi}. ^ ^ 

For the subset {Si} x y the roles of the two fixed points are reversed: {Si} x {Qi} acts as a-limit, 
{Si} X {Ti} as w-limit. 

The fixed point {Pi} x {Qi} is a saddle. While it acts as the a;-limit for the one-parameter family 
of orbits on {Pi} x y, there exists exactly one orbit emanating from it: The ^LRSi orbit Si = 
on E X {Qi} connects {Pi} x {Qi} with {Si} x {Qi}. Analogously, the ""LRSi orbit Si = on 
S X {71} is the unique orbit connecting {Si} x {Ti} with {Pi} x {Ti}. 

Therefore, we have found a two-parameter family of heteroclinic cycles involving the fixed points 
{Si/Pi} X {Ti/Qi}, which can be written down schematically as follows: 

{Si}x{Qi} {Pi}x{Qi} 



infinitely many 



infinitely many (^^) 



{Si} X {Ti} — . {Pi} X {Ti} 

unique 

Evidently, this family of heteroclinic cycles represent a generalization of the anti-LRS cycle (|4T|) . In 
fact, as noted in ((55|) . the flow on {Si /Pi} x y permits alternative completions of the heteroclinic 
cycles: Possible "detours" include additional flxed points, see Figure [9l 

{Si} X {%/T,} {Si} X {Qi} {Pi} X {Qi} {Pi} X {%/T,} 



detour 



I infinitely|many infinitely J many detour (57) 



{Si} X {Q2/Q3} > {Si} X {Ti} > {Pi} X {Ti} > {Pi} X {Q2/Q3} 

detour unique detour 



Sequences of heteroclinic orbits 

As we have seen in Subsection 17. 1[ the flow on the subsets S x {Qi} and S x {71} is represented 
by a family of straight lines that possess a common focal point; see Figure El The anti-LRS line 
Si = is the central line. Orbits that are close to Si = (i.e., orbits where Si is small) connect 
the same sectors as the central line, namely sector (312) with sector (213). Let P(3i2) denote any 
point on sector (312) of dS (the anti-LRS point Pi being an example) and likewise P(2i3) ^-ny 
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{Si} X {Q2} {Si} X {%} 




Figure 9: A depiction of the heteroclinic cycles (|57p: the dashed hnes form a heterochnic cycle of 
the family ((56)) . 



point on sector (312) (such as Si). Symbolically we then write for any orbit on S x {Qi} that is 
close to El =0 



(312) 



{Qi} 



(213) 



{Qi} 



(58) 



and likewise for any orbit on S x {Ti}. For these orbits it is immediate that jSij is monotonically 
decreasing in t. Accordingly, when we invert the direction of time, jSij is monotonically increasing 
along each orbit, so that all orbits diverge from the central line Si — 0. 

These orbits on i7 x {Qi} and S x {71} can be used to construct sequences of heteroclinic orbits. 
Such a sequence of heteroclinic orbits can be written schematically as follows: 



{S{2i3>} X {'^3/2} ^ 



detour 



detour 



{S(213)} X {Ql} 
infinitely 



unique 



many 



{P(312)} X {Ql} 
infinitely 



detour 



{P(312)} X {T3/2} 



many 



detour 



{S(213)} X {Q2/3} 



detour 



{S(213)} X {Ti} 



unique 



{P(312)} X {Ti} 



detour 



{P{312>} X {Q2/3} 

(59) 

Naturally, in the course of the sequence, the points S(2i3) and P(3i2) change; in particular, jSij 
increases (with — r). An example of such behavior is depicted in Figure flOl 

Since jEij grows, any given sequence of the type ([55]) will eventually leave the sectors (213) and 
(312). There exist uncountably many possible continuations; however, while (|59p is associated with 
a unique sequence of Kasner states, the possible continuations of a given sequence are not: They 
give rise to different sequences of Kasner states. Since this makes the analysis difhcult, we refrain 
from going into details. 
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Figure 10: Example of heteroclinic sequence (j59p diverging (as r — > —oo) from a heteroclinic cycle 
of the family ([56)) (represented by a dashed line) . 

7.3 a-limits 

It is obvious that the complexity of the structures that are present on the boundary of the state 
space impedes a rigorous mathematical derivation of the past attractor of orbits in X. We confine 
ourselves to presenting some basic results. 

Theorem [5] implies that the a-limit of every orbit in X must be a subset of dX. Since all fixed 
points on dX act as saddles, there do not exist orbits in X whose a-limit set consists of merely a 
point. In particular, there do not exist solutions that converge to a Kasner solution as t — + 0. 

The past asymptotics of LRS solutions and anti-LRS solutions has been analyzed in Sections [5] 
and [6] The asymptotic behavior of solutions is characterized by oscillations between two different 
Kasner solutions. For solutions that are neither LRS nor anti-LRS we have the following theorem: 

Theorem 5. The a-limit set of a (non-LRS, non-°'LRS ) solution in X comprises a large (probably 
infinite) set of Kasner points. 

Interpretation of the theorem. The theorem states that, as r — > — oo, the solution undergoes an 
infinite sequence of phases ( "epochs" ) , in each of which the behavior of the solution is approximately 
described by a certain Kasner solution. We conjecture that the number of Kasner states the 
(generic) solution passes through in this way is infinite. 

Proof. Consider an orbit 7 in A" which is neither LRS nor ^LRS. Suppose that a point V G dX is 
an element of a(7). Since V G 01(7), the orbit through V and its a-limit point V- and its tj-limit 
point V+ must also lie in a(7). (As a matter of course, since all orbits on dX are heteroclinic orbits, 
V± are fixed points.) Since V+ € ^(t), there exists an orbit emanating from V+ that is contained 
in 01(7) as well; likewise, there exists an orbit converging to V- that lies in 01(7). Continuing in 
this manner we can construct a sequence of fixed points and thus a sequence of associated Kasner 
states ICm n S Z, that is contained in Q!(7). Taking account of the previous analysis of the flow 
on dX it is not difficult to convince oneself that the sequence /C„ does not exhibit any simple 
recurrence (and probably no recurrence at all in the generic case); this is, however, provided that 
V ^ LRSi and V ^ ''LRSi {i = 1,2,3). In that special case, the sequence /C„ is the alternating 
sequence of Kasner states described by (|33p or (|4T]) , respectively. In order to prove the theorem it 
thus remains to show that the intersection of 0(7) with LRSi and ''LRS^ is empty, or, equivalently, 
that 7 cannot converge to any of the heteroclinic cycles ([33]) or (I41|) as r ^ 00. To see that we 
merely note that these heteroclinic cycles are not stable. Consider, for instance, the LRS orbit in 
Cyljjj-j.] . A small perturbation of that orbit results in different a- and w-limit (that are located 
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y3 



yi 




y2 



Figure 11: An orbit converging to a heteroclinic sequence of the type is increasingly forced 
to the boundary of the hexagon by (pT|) . The dashed hues represent the curves 2/2/(1 — 2/2) = 
const 2/3/(1 — 2/3), where the constant is either small (upper dashed line) or large (lower dashed 
line). 



on the Kasner circles, generically) , and thus leads to a quickly increasing deviation from LRS. 
Analogously, consider the anti-LRS orbit = on i7 x {Qi} or S x {%}. Since |Ei| grows along 
the other orbits on i7 x {Qi/Tl}, perturbations of the anti-LRS orbit increase with increasing (— r). 
This establishes the claimed instability of the LRS and ^LRS cycles. □ 



Having shown in the proof of the theorem that there do not exist orbits that converge to the 
LRS/^LRS heteroclinic cycles except for the LRS/^LRS solutions themselves, one follow-up ques- 
tion suggests itself: Do there exist orbits that converge to a sequence of heteroclinic orbits of the 
type (|59p described in the previous subsection? To analyze this question, consider the heteroclinic 
sequence ()59p and suppose that there exists an orbit in X that shadows this sequences in the 
asymptotic regime f = (— r) 00. Along the orbit we thus have 



and accordingly 



d 
df 

2/2 



log 



2/2/(1 - yz) 
2/3/(1 - ya) 



oc 



2/3 



■ exp 



= -6E1 



6 / Eidf 



(60) 



(61) 



i - 2/2 i - 2/3 

Since Ei has a sign along the sequence of heteroclinic orbits, jEidf is monotone in f; in fact, 
since the orbit stays a large amount of time close to the fixed points, the (absolute value of) the 
integral is very large. Consequently, (pT|) implies that 2/2 and 2/3 cannot remain of the same order 
as f increases. Therefore, in the asymptotic regime r — > 00, the orbit is forced to the boundary of 
y, see Figure [TT] In the terminology of the previous subsection we can say that the orbit is forced 
on the detour paths of ([59]). 



A detailed analysis (which we omit here) of the fixed points {P(3i2>} x {Q2/3} and {S(2i3)} x {T3/2} 
reveals that additional branchings can occur at these points. If a > 1 — 2/^/3, then {P(3i2)} x {Q2/3} 
acts as the w-limit not only for one orbit (the orbit {P(3i2)} x {%_} — * {P{3i2)} x {22/3}); but 
for a one-parameter family of orbits. (In fact, {P{3i2)} x {Q2/3} is a sink in Cyl[23i] or Cyl[32i], 
respectively.) It is unclear whether this branching prohibits orbits to shadow sequences of the 
type ISB- 
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The case a < 1 — 2/-\/3 is simpler. If |Ei| is sufficiently small, no branching occurs at the points 
{P{312)} X {Q2/3} and {S(2i3)} X {%,/2]- We conclude that there exist orbits in X that shadow 
the sequence 

{S<213)} X {%/^} {S<213)} X {Ql} {P(312)} X {Q^} {P(312)} X 

detour 1^ J detour 

{S(213)} X {Q2/3} ' {S(213)} X {Ti} ; > {P(312)} X {Ti} > {P(312)} X {Q2/3} 

detour unique detour 

(62) 

for some time in the asymptotic regime r — > — cx). Accordingly, if becomes sufficiently small 
in the asymptotic evolution of a solution, then the solution enters an "anti-LRS phase" (where its 
behavior resembles the behavior of anti-LRS solutions). The length of this phase (as measured, 
e.g., by the number of oscillations in the sequence) is inversely proportional to the initial value of 
I Si I and can thus be arbitrarily long. During an anti-LRS phase, jSij increases until either the 
orbit branches off at {P(3i2)} x {Q2/3} or {8(213)} x {'?3/2} when |Si| has become large enough, 
or the orbit leaves the sectors (213) or (312). 

It remains to ask whether more general considerations than (j61l) can lead to a further exclusion of 
parts of dX as possible a-limit sets. For instance, it might turn out to be true that the a-limit of 
a generic orbit in A" is a subset of the set 

U = (^^ X aj;) u y (r X {qj) u |J (r x {%} 

ijk i i 
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